Dynamic stability of spindles controlled by molecular motor kinetics 
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We analyze the role of the force-dependent kinetics of motor proteins in the stability of antiparallel 
arrays of polar filaments, such as those in the mitotic spindle. We determine the possible stable 
structures and show that there exists an instability associated to the collective behavior of motors 
that leads to the collapse of the spindle. Our analysis provides a general framework to understand 
several experimental observations in eukaryotic cell division. 
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Living cells display several structures that arise from 
the self-organization of polar filaments and motor pro- 
teins Q. Several in vitro studies have shown that 
mixtures of kincsin motors and microtubules (MTs) 
can spontaneously develop complex spatio-temporal pat- 
terns 2$. These self-organization processes are essential 
for eukaryotic cell division jjj]. During mitosis, motor 
proteins organize MTs in a bipolar structure, the mitotic 
spindle, which serves as a scaffold to transmit the neces- 
sary forces for chromosome segregation 0. The spindle 
consists of two MT asters that overlap in the central re- 
gion. The MTs, with their minus-ends located at the 
aster poles, are crosslinked by many different motor pro- 
teins 0,0- One particular type of motors, the plus-ended 
bipolar kinesins (e.g. Eg5 or Klp61F), has been shown 
to be essential for the spindle stability. A decrease in 
their concentration below a certain threshold causes the 
spindle collapse 0, @ , and their total inhibition prevents 
bipolar spindle formation ||. In addition, Eg5 motors 
have been shown to drive the MT poleward flux |7( and 
homolog motors to induce the formation of (interpolar) 
MT bundles @. 

Bipolar motors are composed of two connected units, 
each one composed of two motor domains. Both units 
can move simultaneously and independently on MTs • 
These motors are able to crosslink MTs and slide them 
with respect to each other when they are in an antipar- 
allel configuration like in the central region of the 
spindle (Fig. |JJi,c). As a result, these motors produce an 
outward force along the spindle axis and generate a MT 
flux toward the poles 0. Typical forces involved in mi- 
tosis lay in the nanoNewton range . Since individual 
motors cannot exert forces larger than a few picoNew- 
tons, their collective action is required to ensure the sta- 
bility of the mitotic spindle. At metaphase, this dynamic 



structure reaches a steady state with MTs of nearly con- 
stant length undergoing permanent treadmilling 0, Il2| , 
end and depolymerizing at the — 



polymerizing at the 
end. 

The theoretical study of motors and MT mixtures has 
been recently addressed using continuum coarse grained 
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FIG. 1: (a) Antiparallel array of MTs under the action of 
a longitudinal force F. Buckled and non-buckled arrays are 
shown. The minus and plus ends of MTs are depicted as — 
and + respectively. Motors are represented by dots. In the 
non-overlapping region of length L there are no antiparallel 
filaments and the motors are not subject to any force. The 
motors in the overlapping region of length I sustain the struc- 
ture by crosslinking and sliding antiparallel filaments, (b-c) 
Possible kinetic events of a motor in the non-overlapping re- 
gion (b) and in the overlapping region (c). The velocities V 
and V axe the crosslinking motor velocity and the MT pole- 
ward velocity respectively. 



descriptions [l3L IT^ . which have elucidated their basic 
self-organizing principles. However, the coupling between 
force-dependent motor kinetics and local forces in self- 
organized structures has not been addressed. In this Let- 
ter, we study the dynamic stability of antiparallel arrays 
of MTs under the action of longitudinal forces, in the 
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presence of molecular motors able to collectively hold the 
structure by stochastically crosslinking the filaments. We 
analyze the effects of the motor kinetics on the stability 
of the structure, and show that several phenomena ob- 
served in eukaryotic cell division can be understood on 
the same physical basis. 

In order to comprehend the basic physical mechanisms 
controlling the stability of a spindle, we concentrate on 
a simplified geometry. We consider a pair of antiparal- 
lel MTs (or an antiparallel MT bundle) of fixed length, 
under the action of an inward force F (Fig. ^,). In 
the steady state, there is a region of length I where the 
antiparallel filament array overlaps (overlapping region) . 
The motors in this region can crosslink antiparallel MTs 
and slide them in opposite directions, generating an out- 
ward force that balances the applied force F. We assume 
the antiparallel MT sliding to be the only mechanism 
generating the poleward MT flux, as suggested experi- 
mentally [3 . Out of this overlapping region there are two 
regions of length L (non-overlapping regions; see Fig. ^i) 
where motors cannot apply forces to sustain the spindle. 
Since we concentrate on the spindle stability and do not 
address the mechanisms that determine L and I, we take 
them as given parameters. 

The motors in the non-overlapping region can be either 
bound to a MT or freely diffusing in the bulk (Fig. ^d). 
We assume a constant bulk motor density, p 3 D [la], ar± d 
consider the motors in the bulk to attach onto MTs at a 
rate k% D ■ Once bound to a MT, a motor moves convec- 
tively with a mean velocity Vq toward the plus-end of the 
MT and detaches at a rate At mean field level, the 
dynamics of the bound motor density, p(s,t), can then 
be expressed as 0, 0] 

d t p(s, t) + d s J(s, t) = -k° u p{s, t) + kl D ■ 
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where s is the position along the MT as measured from 
the MT minus ends and J(s,t) is the flux of bound mo- 
tors. For simplicity, we assume the bound motors to be in 
a low density phase and write J(s, t) = p(s, t)(Vo — V) in 
the laboratory reference frame, with V being the velocity 
of MTs toward the poles (Fig.^D,c). 

In the central overlapping region, bipolar motors can 
be either in a crosslinking state or in a bound state. 
In the former state both motor units are attached to 
a pair of antiparallel MTs, sliding them in opposite di- 
rections and supporting a fraction of the total force F. 
As a result, the crosslinking motors move with a force- 
dependent velocity V (Fig. QJ;). Based on experimen- 
tal observations , we write a linear force- velocity re- 
lation, V — Vb (1 — fm/fs), where f m is the load ap- 
plied on the motor and f s its stall force. We consider a 
number n c of independent crosslinking motors to equally 
share the total applied force, so that f m = F/n c . As 
the poleward MT movement is driven only by these mo- 
tors, we identify V = V (Fig. [TJ;). Each unit of a motor 
in the crosslinking state can detach at a force-dependent 



rate k u (f m ) = k°exp(f m b/K B T) (Kramers theory 0), 
where b is a length in the nanometer scale characteriz- 
ing the activated process and KbT the thermal energy. 
Such exponential sensitivity to applied load has indeed 
been observed experimentally |20j . After the detachment 
of one motor unit, the bipolar motor is only bound to one 
MT and unable to apply force. Such motor can either de- 
tach the bound motor unit left at a rate and diffuse 
into the bulk, or re-attach the unbound motor unit at 
a rate k\, and become a crosslinking motor again. The 
motors in the bulk can also attach directly to the MTs 
in the overlapping region at a rate k\ D . 

The relevant variables being the number of motors sus- 
taining the spindle, we neglect their spatial distribution 
in the overlapping region. Accordingly, the equations for 
the average number of crosslinking and bound motors, 
n c and rib respectively, read 

^ = k b n b - k u (n c )n c , (2) 
at 

drib 
~~dt 



= 2J(L, t) + kl D p 3D £ + k u (n c )n c - (k° + k b ) n b , 



where J{L, t) is the convective flux of bound motors com- 
ing from a non-overlapping region, with L being the ar- 
clength of a MT from its pole to the overlapping region 
(Fig. nji). The value of J(L,t) is determined by the so- 
lution of Eq. ^ When the motor processivity length, 
lp = Vo/k-u, is smaller than the characteristic spindle 
length (l p <C L [2l|), the flux J(L,t) is determined by a 
constant bound motor density, p = k^ D p^^/k^, fixed by 
the exchange of motors with the bulk. In this case the 
dynamics of n c and n b are decoupled from p close to the 
pole. 

The existence of antiparallel MT arrays under an ex- 
ternal load F is determined by the balance between mo- 
tor attachment and detachment fluxes, as given by the 
steady state solutions, {nl,^}, of Eqs.^13 which read 
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The quantities h c = n c /F and hb = rib/F are the nor- 
malized numbers of crosslinking and bound motors re- 
spectively. The ratio, F = Fb/KsT, between the force 
F and the characteristic detachment force KsT/b sets 
the natural scale of motors. The relevant dimensionless 
lengths are L = aL/l p and I = a£/l p , and the param- 
eter a = fsb/KsT quantifies the sensitivity of motor 
detachment to force. The asymmetry in motor attach- 
ment/detachment events at vanishing load is character- 
ized by 7 = h/k°. 

There always exists a critical value, S m , below which 
there are no solutions of Eq. [3] This situation corre- 
sponds to an attachment flux of crosslinking motors that 
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FIG. 2: (a-b) Dependence of the critical value S c on L and I in 
the limiting cases where (a) < < 1 and (b) I > 1 respectively, 
(c) Possible spindle structures as the bulk motor density p^o 
and the force F are varied {psn = S c a(k l l) 3 / k^ D Vokb sets 
the density units). Above p^o n {F), buckled (straight) stable 
spindles exists for F > F B (F < F B ). Below p™£ n (F) no 
stable spindles exist. 



can not balance their detachment flux, leading to the 
loss of all crosslinking motors and inducing the spindle 
collapse. Associated to the critical point S m , there is a 
minimum number of crosslinking motors, n™, necessary 
to sustain a spindle, whose value is given implicitly by 
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The actual value of S m is obtained by substituting n' c in 
Eq.GDby the solution of Eq.H 

In order to determine the stability of the structures, 
we perform a linear stability analysis of the solutions of 
Eq. For L ^> l p it can be shown that the fluctuations 
in p(s) are negligible and the spindle stability depends 
only on the dynamics of n c and nt, [13] • Stable spindles 
exist above a critical value S c . The precise expression for 
this critical point S c depends on the ratio between mo- 
tor attachment/detachment rates at vanishing load. For 
7 > 7 C = exp (1/n™) (1 - n™ 1 ) /n" 1 - 1, the transition 
from an unstable array to a stable spindle corresponds 
to a saddle-node bifurcation at 6 C = S m . On the other 
hand, if 7 < j c , this transition corresponds to a global 
bifurcation (saddle-connection J23|) at a value S c > S m . 
Regardless the value of 7, the same qualitative scenario is 
observed as S is varied and we restrict the following dis- 
cussion to the regime 7 > j c without loss of generality. 
The threshold value, 5 Cl is fixed only by L and £. Typi- 
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FIG. 3: Steady state analytical solutions for (a) the number 
of crosslinking motors and (b) the velocity of the crosslinking 
motors (sliding velocity), as a function of the bulk density of 
motors. The different curves represent different values of I in 
the limiting regime I > 1. 



cally L ^> l p and S c becomes independent of L (Fig.^Jj), 
while I modifies 5 C slightly (Fig. ^jp). In this limit, the 
spindle morphology and its stability are decoupled. As 
the motor properties and the lengths L and £ are difficult 
to modify experimentally, S appears as the natural con- 
trol parameter for the spindle stability, since it depends 
both on the applied force, F, and the bulk motor density, 

P3D- 

The existence of a critical value S c implies that for a 
MT array under the action of a load, F, there exists a 
minimal motor bulk concentration, p™p n , below which 
no stable spindles are found. Using the definition of 5 
(Eq. |2| , this minimal motor density reads 



min _ u " u c ^ 
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When p 3 £> > the stable spindle may be either 

straight or buckled depending on the value of the com- 
pressive force, F. For forces below (above) the buckling 
force Fb = (it / 2) 2 B / L 2 of the structure (B being its 
bending rigidity), the stable spindle is straight (buck- 
led). In Fig.[2t we plot the structures that can be found 
as the bulk density of motors and the force applied on 
the spindle are varied. Indeed, recent experimental ob- 
servations have shown that the progressive inhibition of 
Eg5 motors leads to the collapse of the spindle at a finite 
bulk motor density |3- Moreover, the total inhibition 
of homolog motors (Klp61F) has been shown to prevent 
bipolar spindle formation in vivo [8| , in accordance with 
our predictions. 

Above the threshold density p™^ 71 , there is a finite 
amount of crosslinking motors n c collectively holding the 
spindle. In Fig. [3^ we plot the stable solutions of Eq.|31as 
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a function of the bulk concentration of motors. Increas- 
ing values of p-^D and £ leads to larger motor attachment 
fluxes, that result in a larger amount of crosslinking mo- 
tors. For a living cell in native conditions, the MTs in the 
mitotic spindle are typically buckled 0. Therefore, the 
force applied on the MTs is of order Fb which, for single 
MTs (5/im in length) is about lpN. Using this value for 
the force F, the number of crosslinking motors leading 
to a stable antiparallel array turns out to be very small 
(~ 2). In this case fluctuations would dominate and, al- 
though stable arrays could be transiently formed, their 
lifetime would be too short (on the time scale of motor 
detachment). Since the buckling force of a MT bundle 
can be at least one order of magnitude larger [24| , stable 
interpolar MT bundles require tenths of crosslinking mo- 
tors and provide robust spindles with lifetimes over the 
time scale of the division process. Interpolar MT bun- 
dles are indeed observed in several organisms during cell 
division 0, EE E3| - 

The speed of the MT flux toward the poles is deter- 
mined by the MT sliding velocity, V, given by the velocity 
of the crosslinking motors. In Fig. [3t> we represent this 
sliding velocity as a function of p^o, for different values 
of £. It decreases from its maximal value Vq as the bulk 
motor density is decreased, and it is typically finite for 
the minimal density p™™ at which the spindle collapses, 
as observed experimentally • At high motor concentra- 
tions psd 3> the crosslinking motors move nearly 
at their maximal velocity Vo and the MTs move poleward 
at this velocity consequently. As the motors in the non- 
overlapping regions move at velocity Vo with respect to 
the MTs in the spindle (Fig.[TJ)), they appear static in the 
laboratory reference frame. This apparent motor stillness 
has indeed been observed experimentally Our anal- 
ysis predicts that decreasing the bulk motor density p^o 
would allow the observation of motor movement in the 
spindle. This observation would provide further insight 
on the understanding of the mitotic spindle structure. 

The present approach highlights the importance of 
force-dependent motor kinetics on the self-organization 
of MTs and motors. In particular, we show that mo- 
tor kinetics is a key factor in the stability of spindle-like 
structures under applied load. Moreover, we have de- 
veloped a framework through which several observations 
in eukaryotic cell division, namely the spindle collapse, 
the MT poleward flux, the static appearance of motors 
in the spindle, and the existence of MT bundles, can be 
explained on a common physical ground. The observa- 
tion of motor motion upon a decrease of motor density, 
as predicted here, would provide strong evidence for the 
self-organization of motors and MTs as the underlying 
principle of mitotic spindle assembly. 
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